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Abstract 

Radiative transfer modelling of AGB circumstellar envelopes is applied 
to a sample of AGB stars previously observed with the mid-IR imaging 
camera TIRCAM (Busso et al. 1996: Paper I). We present the results of 
our simulations, aimed at deriving the physical parameters of the envelope, 
such as the optical depth and the radial thermal structure, the mass loss 
and the dust-to-gas mass ratio. The chemical composition of the dust in 
the observed envelopes is discussed. The ability of different sets of dust 
opacities to fit the mid-infrared spectra is evaluated. The hypothesis of 
dust grain aging and annealing in O-rich envelopes is considered in order 
to explain an apparent inadequacy of the availabe opacities to describe 
the variety of observed spectra, as previously noted by other authors. 
Various possible origins of the discrepancies are discussed, together with 
their consequences on the dust grain formation processes. Keywords: 
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1 Introduction 



Asymptotic Giant Branch (AGB) stars are characterized by intense mass loss 
processes that play a central role in the star's subsequent evolution and lead to 
the formation of a cold circumstellar envelope (CSE) of gas and dust. The enve- 
lope chemical composition is determined by the amount of carbon enrichment in 
the stellar atmosphere, due to mixing processes (the so called third dredge-up), 
induced by thermal pulses near the end of the AGB phase. C-rich envelopes 
(having [C]/[0]>1 by number) mainly exist around carbon stars, while O-rich 
CSEs ([C]/[0]<1) are associated to M giants. The transition from O-rich to 
C-rich envelopes is controlled by the extent of the third dredge-up, and it can 
occur only for AGB stars in a narrow mass range: indeed, stars with a core mass 
lower than 0.6 M (Straniero et al. 1995) will not undergo the third dredge-up, 
while H burning at the base of the convective envelope (Hot Bottom Burning) 
is expected to deplete carbon in the envelope of stars more massive than ~ 
5 Mq (Wood et al. 1983; Boothroyd et al. 1993; Frost & Lattanzio 1995), 
preventing them from becoming C-rich. 

Despite the low dust-to-gas mass ratio (in the range ~0. 001-0. 01), the opti- 
cal properties of AGB CSEs are mainly determined by the dust grains: silicates 
in O-rich envelopes (Pegourie & Papoular 1985) and a mixture of Hydrogenated 
Amorphous Carbon (HAC, Jones et al. 1990; Duley 1993) with inclusions of 
SiC (Skinner & Whitmore 1988) and possibly Policyclic Aromatic Hydrocar- 
bons (PAH, Puget & Leger, 1989; Cherchneff & Barker 1992) in C-rich CSEs. 
Both oxidic and carbonaceous star dust is characterized by vibrational bands 
positioned in the mid-IR window (at 9.7 and 18 /xm for silicates, at 11.3 (im for 
SiC and at 3.3, 6.2, 7.7, 8.6 and 11.3 /im for PAH); furthermore, most of the 
continous thermal radiation from the dusty optically opaque envelope is emitted 
in the same wavelenght range. 

For these reasons a mid-IR search was carried out with the imaging camera 
TIRCAM (Busso et al. 1996, hereafter Paper I) for a sample of 16 sources, 
including both AGB (O-rich and C-rich) and post-AGB stars. In Paper I, 
photometry and colors were derived using 10% bandwidth filters at 8.8, 9.8, 
11.7 and 12.5 fim, and compared with the Low Resolution Spectra (LRS) from 
IRAS (1986), in order to establish suitable photometric criteria for discrimi- 
nating between O-rich and C-rich sources and for estimating mass loss rates. 
The observations formed the basis for discussing the evolutionary status of the 
sources and their mass loss history. In Section 2 we discuss the general features 
of their spectra, from the point of view of the optical properties that can be 
inferred for dust. 

Detailed numerical simulations of the CSE thermal structure is the subse- 
quent step for a better determination of the physical parameters. One has to 
compute the source spectrum and the radial brightness distribution for a wide 
range of input parameters; comparison with observations will then allow us to 
estimate the optical depth of the envelopes, to derive mass loss rates, and to 
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extract information on chemical abundances which are relevant to study the 
stellar nucleosynthesis (see e.g. Busso et al. 1995). With these aims in mind, 
we developed a numerical code for solving the problem of non-grey radiative 
transfer through a spherically symmetric dust shell around an evolved AGB 
star. We present our model in Section 3, with emphasis on the input model 
parameters and on the dust opacities used. In Section 4 the model results are 
shown for the sources observed in Paper I whose LRS is known, and a new 
estimate for the dust mass loss rates and the dust-to-gas mass ratio is obtained. 

In Section 5 a more general discussion of the model spectra is carried on. 
We point out how available dust opacities seem not to be completely adequate 
for a detailed simulation of the mid-IR features of circumstellar dust, as pre- 
viously suggested by different authors (Little-Marenin & Little 1990, hereafter 
LML90; Simpson 1991), and we test possible explanations. Our conclusions are 
summarized in Section 6. 

2 Spectra of AGB Circumstellar Envelopes 

The IRAS color-color diagram (van der Veen & Habing 1988) obtained by IR 
photometry at 12, 25 and 60 «m, provides a first-order discrimination of AGB 
sources, according to their infrared excess, variability and C abundance. How- 
ever, as shown in Paper I, such IRAS color-color diagram can lead to misclas- 
sification of AGB CSEs due to the large width of IRAS filters (this is true in 
particular for the one at 12 am, which cannot discriminate between silicates and 
carbonaceous features). Although narrow band filters allow a better discrimi- 
nation (as in the case of TIRCAM photometric system), a detailed description 
of the dust emission/ absorption features requires mid-IR spectra. 

IRAS LRS (1986) provide spectral energy distributions in the wavelength 
range 7.7-22.6 /zm for about 4,000 AGB sources, which are classified with a 
two-digit code. The first digit (main class) is related to the spectral index [3 
of the continuum (assuming F\ oc X 13 ), while the second indicates the strength 
of the main feature, in emission or absorption, present in the mid-IR window. 
In particular, spectra of O-rich sources are identified by the 9.7 /im silicate 
feature, that is found in emission or absorption in the LRS classes 2n and 3n, 
respectively. C-rich envelopes are mainly found in the An class, defined by the 
presence of the 11.3 um SiC feature in emission. Post-AGB objects of both 
types are in the main classes 6n and 7n, characterized by very large infrared 
excess. Finally, both O-rich and C-rich sources may fall in classes In and On 
(featureless spectra), where a few LRS with self-absorbed SiC feature can be 
found (Omont ct al. 1993). 

Table 1 lists the spectral characteristics of the sources in Paper I with known 
LRS; the star AFGL 1822 was added to the TIRCAM sample in order to have 
an evolved O-rich envelope with the 9.7 um silicate feature in absorption. Post- 
AGB sources have not been included because their modelling involves very dif- 
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II 

# 


iKAb name 


Name 


Type 


ar type 


T no 


Lixo sub. 


T 1 


02143+4404 


W And 


S-star 


S8,2c(M) 


22 


S 


2 


08525+1725 


X CnC 


C-rich 


C5,4(SRb) 


42 




3 


09425+3444 


R LMi 


O-rich 


M7e(M) 


24 


Sil+ 


4 


09452+1330 


CW Leo 


C-rich 


Ce(M) 


43 




5 


12427+4542 


Y CVn 


C-rich 


C5,5J(SRb) 


42 




6 


13001+0527 


RT Vir 


O-rich 


M8III(SRb) 


21 




7 


15193+3132 


S CrB 


O-rich 


M6,5e(M) 


24 


Sil+ 


8 


15255+1944 


WX Ser 


O-rich 


M8e(M) 


29 


Sil 


9 


16029-3041 


AFGL 1822 


O-rich 


M 


32 




10 


16235+1900 


U Her 


O-rich 


M7IIIe(M) 


23 


Sil+ 


11 


18040-0941 


FX Ser 


C-rich 


C(Lb) 


44 




12 


18397+1738 


NSV 11225 


C-rich 


Ce 


43 





Table 1: List of AGB sources from Paper I with known LRS. The evolved O-rich 
envelope AFGL 1822 was added to have a in LRS class source. Last column: 
LML90 LRS subclass, when available. 

ferent physical conditions, dominated by the hot radiation emitted by the central 
star (B or A spectral type) on its way to become a white dwarf. 

Note that the LRS class 2n contains sources with different characteristics; a 
tentative classification is in the LML90 paper, where the silicate band profiles 
in optically thin O-rich envelopes are divided into seven subclasses, according 
to their shape and the presence of secondary peaks. The subclass for our O- 
rich sources, when given by LML90, is reported in Table 1; the S-star W And 
is tentatively attributed to the "S" feature subclass, characterized by a broad 
peak at 10.3 fim. The other sources fall in subclass "Sil" (with a strong feature 
peaked at 9.8 /*m, typical of 2n LRS with n > 4) and subclass "Sil+" (having a 
secondary "bump" at 11.3 /im). LML90, and Simpson (1991) suggest that these 
differences be related to a kind of "mincralogical diversity" in the oxidic CSE 
dust, although this explanation is still controversial (see e.g. Ivezic & Elitzur 
1995, hereafter IE95). We tried to model the shape of the silicate feature in our 
sources belonging to different subclasses using several sets of opacities and dust 
mixtures. 

3 Our Model 

Radiative transfer in circumstellar dust shells has been previously discussed by 
many authors (Rowan-Robinson 1980; Griffin 1990; Justtanont & Tielens 1992, 
Groenewegen, 1993 and, more recently, Hashimoto 1995; a detailed discussion 
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of these models is given in Habing, 1996). A fully self-consistent modelling of 
envelopes around evolved late-type stars requires coupling of radiative transfer 
with the hydrodynamic equations of motion for the two interacting fluids (gas 
and dust) which constitute the stellar wind. Such an approach has been at- 
tempted by IE95: their treatment, in the hypothesis of a steady-state outflow 
due to constant dust-driven stellar wind, can successfully account for the ob- 
served mass loss rates and is able to reproduce the IRAS colors and visibility 
functions (see Ivezic and Elitzur 1996a, b) for a large sample of AGB objects in 
all LRS classes. More recently, models taking into account time dependent out- 
flows have been developed to study the effect of sudden variations in the stellar 
parameters in connection with thermal pulses (Shonberner 1996). This analysis 
is necessary in order to follow the latest stages of AGB evolution, which lead to 
the formation of a planetary nebula; however, time variations of the mass loss 
rate, on a thermal pulse time scale, have negligible effects on the IR properties 
of at least 95% of late-type stars associated with IRAS sources (IE95). 

In our model we totally neglect the dynamics of the envelope; rather, we 
concentrate our efforts on obtaining an accurate fit of the observed sources, 
paying special attention to the mid-IR interval, where the TIRCAM spectral 
bands are located (see Paper I). The model parameters are fixed by fitting 
the IRAS and TIRCAM data. For any dust mixture, the spectrum (and in 
particular the shape of the features) is used to estimate the optical depth t„ of 
the envelope, proportional to the dust opacity and the mass loss rate. Near- and 
Mid-IR photometric data allow us to determine the stellar parameters T e f f , i?„ 
(the central star is assumed to radiate a black body spectrum at temperature 
T e //), and the inner CSE radius J?i; the outer radius R2 is derived by fitting 
the IRAS far-IR photometry at 60 and 100 /zm. The source distance enters in 
the model only as a scale factor for the flux densities, and is determined by 
normalizing the spectra with the observed photometry. We first assumed the 
value estimated by Loup et al. (1993); then we computed our own value of 
d = y/L,lF lkpc , where L* = 4naRlT* ff is the total luminosity and Fik pc the 
total flux (in L©, assuming a reference distance of 1 kpc) obtained by Loup et 
al. (1993) using IRAS fluxes. 

The code computes iteratively a self-consistent thermal structure of the en- 
velope; the emergent spectra is calculated taking into account the effect of non- 
isotropic scattering, absorption and thermal reemission by grains. The compu- 
tation is performed on the hypotheses of: (i) spherical symmetry of the dust 
shell, with an n(r) oc r~ 2 density distribution, consistent with a steady outflow 
at constant velocity, (ii) balance between absorption and emission by the dust 
grains, and (iii) LTE dust radiation at the local temperature T(r). 

The dust grain composition is simulated adopting different sets of dust opac- 
ities Q v \ the dust grain size is usually assumed to satisfy a size distribution 
n(a) oc a~ 3 - 5 with 0.01 ^m < a < 0.25 ^m (Mathis, Rumpl and Nordsieck 1977, 
hereafter MRN; see e.g. Jura, 1996 for more recent data). The main effect 
of considering a grain size distribution is a broadening of the 9.8 /an silicate 
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feature (Simpson 1991), due to the non linear dependance of the opacity on a: 
this is however important only for A < a and in the case of AGB CSE grains 
can be neglected. Hence we simply adopted for the grain size the average MRN 
value a = 0.1 /im. 

3.1 Dust Opacities 

Following a Chan & Kwok (1990) scaling theorem, the solution of the radiative 
transfer equation in spherical symmetry is scaled on the flux-averaged optical 
depth Tp = J Q T v F v dvj J Q F v dv, where: 



In equation 1, t v depends on frequency only through the opacity Q v , which 
is thus the main parameter for modelling the spectra; for this reason, a great 
part of our preliminary work consisted in testing the various opacity profiles 
available in the literature. We now present the results of our analysis, separately 
for O-rich and C-rich dust. 

It is well known that O-rich dust contains essentially silicate grains: they can 
be present as piroxens (Mg,Fe)SiC>3, or as olivinae ((MgjFe^SiCU). Actually, 
the olivina opacity profile of Kratschmer & Huffmann (1979) is suitable for 
supergiants, but not for Miras like WX Serpentis (see Griffin 1993). In the case 
of circumstellar envelopes the silicate dust probably consist of piroxens only, 
with various types of impurities making them "dirty silicates" . We considered 
five opacity profiles for astronomical silicates. 

The first is the one by Jones & Merril (1976); they introduced the concept of 
"dirty silicates" , though many details of their opacity profile were subsequently 
modified. 

The second is by Draine & Lee (1984); despite its wide use and accuracy, 
this is not well suited for our needs, since it is based on observations of the 
interstellar medium, and not of the circumstellar one (see also Ossenkopf et al. 
1992, hereafter OS92, and Griffin 1993). 

The third is by Volk & Kwok (1988): it is based on the observed spectra 
of 467 AGB stars with good quality photometric data in all four IRAS bands; 
the spectrum obtained from this opacity is similar to the observed one, but the 
9.8 \xm feature is slightly too much asymmetric and its peak is shifted towards 
redder wavelengths. 

The fourth profile is the one by OS92: it is based on the previous but it is 
chemically and physically consistent, in the sense that computes the complex 
dielectric function for the dust using the Kramers-Kronig relations along with 
the Mie theory. Taking into account the possibility of dust annealing in the 
CSE, as illustrated by Stencel et al. (1990) and Nuth and Hecht (1990), the 
authors give the opacity in two versions, one for "warm" silicates (hereafter 
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Ossl) and the other for "cold" ones (hereafter Oss2), depending on the physical 
conditions at their condensation. 

The last and more recent is by David & Pegourie (1995), and is computed 
using Mie theory for a set of more than 300 IRAS sources with known LRS. 

We obtained the best results when fitting our AGB O-rich sources with OS92 
Ossl and Oss2 opacities. 

In the case of C-rich envelopes, the continuum in the spectra is usually 
ascribed to the presence of graphite or amourphous carbon, while the charac- 
teristic band at 11.3 /xm is commonly attributed to SiC, although the presence 
of PAH has been proposed as well. Graphite is now commonly ruled out, be- 
cause of the absence in the observed spectra of its typical narrow band at 11.52 
/im (Draine 1984); moreover, the whole shape of the spectra obtained by us- 
ing graphite (the most used opacities are by Draine & Lee 1984) does not fit 
the observed ones, as we tested directly. Therefore the main constituent of 
CSE carbonaceous dust seems to be amourphous carbon: due to the variabil- 
ity in its physical and optical properties, a comparison with observational data 
is fundamental. An accurate enough opacity profile, observationally tested for 
amorphous carbon is the one by Martin & Rogers (1987). More recently Rouleau 
& Martin (1991) presented a detailed study giving accurate sets of optical con- 
stants derived by Bussolctti et al. (1987): the set named AC1, in particular, is 
similar to the Martin & Rogers (1987) opacity, although more detailed; this is 
the opacity we used for our simulations. 

Finally, in order to obtain spectra with the characteristic 11.3 ^m band, we 
admitted the presence of SiC mixed with the amorphous carbon: we choose the 
opacity set for a-SiC (the exagonal rhombocdric cristalline form of SiC) from 
Pegourie (1988), which at the moment is the only accurate profile available in 
the literature; to obtain the optical properties for the mixture, we simply made 
a weighted average of the two species: we never found necessary to take a SiC 
percentage greater than 10%. This strongly depends on the opacity chosen for 
amorphous carbon; IE95, for example, used a more absorptive kind of amor- 
phous carbon, together with Pegourie (1988) data for SiC, and found a SiC 
percentage up to 30%. 

In Figure 1 we plot the final opacity profiles for silicates, amorphous carbon 
and SiC used in our simulations. 

In our simulations we have not introduced Polycyclic Aromatic Hydrocar- 
bons (PAH), because they are mainly excited by UV radiation, and their con- 
tribution to the dust heating and the final spectra is low for AGB stars, at least 
if a hot star companion is not present (Buss et al. 1991). 

4 Model Results 

In order to fit the IR spectra of the sources in our sample, for each set of dust 
opacities we computed a grid of models parameterized on the 10 /im optical 
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Figure 1: Dust opacities used in our computations: Ossenkopfl "dirty silicates" 
profile by OS92 (dot-dash-dot line), amourphous carbon AC1 by Rouleau & 
Martin 1991 (dashed line) and Pegourie (1988) SiC (dotted line). 



depth tio , which is proportional to the dust mass loss rate of the envelope: 

167T Pdd „ 

M d = —^-r w R 1 v e (2) 

where pd is the density of the dust grains (we assumed 3.0 g/cm 3 for OS92 
silicates, 1.85 g/cm 3 for AC1 amorphous carbon and 2.5 g/cm 3 for Pegourie 
1988 SiC), Qio is the opacity at 10 /<m, and v e is the outflow velocity of the 
dust. The models were computed assuming d = 1 kpc, i?* = 10 13 cm, T* = 
2500 K, i? 2 = 2000 J?i, iii = 4 i?* if n0 < 1 and R x = 5 R« if t w > 1. 

Figure 2 shows model spectra (in the wavelength interval 1-100 /im) of O- 
rich envelopes for increasing logrio, from —3 to 1; in Figure 3 the sequence of 
C-rich models with logrio from —3 to is plotted. In both cases greater optical 
depths cause larger infrared excesses, and the dust features (silicates 9.8 yum 
and SiC 11.3 fim) change from emission to absorption. We stopped the C-rich 
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# 


Name 


tio 


v e 


M 


M d 




L, 




d 








km/s 


M /yr 


M Q /yr 


% 


L© 


L0 


kpc 


1 


W And 


0.17 


20. 


9.7 • 10"'' 


2.5 ■ 10~ a 


0.26 


2.9 • 10 4 


4.4 • 10* 


0.26 


2 


X CnC 


0.15 


9.7 


4.6 • 10"'' 


1.5 ■ 10-° 


0.33 


1.3 • 10 4 


2.1 • 10* 


0.79 


3 


R LMi 


0.15 


7.0 


5.0 • 10 _Y 


1.0 ■ 10~ a 


0.20 


5.7 • 10 3 


9.4 • 10 4 


0.25 


4 


CW Loo 


1.30 


15. 


4.8 • 10~ : > 


5.0 ■ 10~ a 


0.10 


4.6 • 10 4 


6.6 • 10 b 


0.26 


5 


Y CVn 


0.01 


8.2 


2.8 • 10~ Y 


1.7 ■ 10~ lu 


0.06 


6.0 • 10 3 


8.7 • 10 4 


0.26 


6 


RT Vir 


0.15 


9.3 


1.3 • 10~ b 


4.3 ■ 10~ a 


0.33 


5.9 • 10 4 


1.3 • 10 b 


0.67 


7 


S CrB 


0.15 


6.3 


6.0 • 10 _Y 


8.0 ■ 10~ lu 


0.14 


4.6 • 10 3 


2.3 • 10 4 


0.45 


8 


WX Scr 


1.30 


8.8 


5.3 • 10"'' 


8.2 ■ 10~' J 


1.55 


3.2 • 10 3 


7.2 • 10" 


0.67 


9 


AFGL 1822 


8.00 


17. 


2.5 • 10"" 


2.0 ■ 10~ s 


0.08 


1.2 • 10 4 


2.8 • 10 J 


2.10 


10 


U Her 


0.15 


13. 


2.6 • 10~ y 


1.1 ■ 10~ a 


0.42 


1.1 • 10' 1 


3.2 • 10 b 


0.06 


11 


FX Scr 


0.20 


26. 


1.2 • 10~ 5 


1.7 ■ 10~ s 


0.14 


8.1 • 10 3 


8.4 • 10' 1 


0.98 


12 


NSV 11225 


0.02 


14. 


9.1 • 10~ B 


1.1 ■ 10~ s 


0.12 


1.2 • 10 4 


2.8 • 10 4 


0.65 



Table 2: Model estimates for sources in Table 1. Envelope optical depth Tio, 
dust mass-loss Md, dust to gas ratio /z, total luminosity L* and distance d for 
each source are determined by model fitting parameters, while total mass loss 
M, total fluxes Fik pc and gas outflow velocities v e are derived by IR-radio- 
millimetric observations from Loup et al. 1993. 

sequence at no — 3 because this can be taken as a limit for this class of sources, 
considering that SiC features in absorption are rare in the LRS database (see 
Omont ct al. 1993). 

The van der Venn & Habing (1988) IRAS color-color diagram for the com- 
puted model sequences is shown in Fig. 4. The two curves are associated to 
increasing optical depth, and thus to increasing mass loss; note that they cross 
the appropriate regions of the diagram for the two classes of sources. The curves 
should not be considered as evolutionary tracks, because there is indication that 
mass loss rates do not increase monotonically during AGB evolution, rather 
they can have sudden variations in connection with thermal pulses (Vassiliadis 
& Wood 1993; see also Paper I). 

4.1 Spectral Distribution 

By means of our model we finally fitted the spectral distribution of the sources 
listed in Table 1 and we obtained some relevant physical parameters for the star 
(T e ff, i?*) and the envelope (i?i, i?2, tio)- 

Moreover, we derived the dust mass loss Md (using Eq. 1, with v e from Loup 
et al. 1993) and the bolometric luminosity L to t — 4tt<tR^T^ j. 

From Loup et al. (1993) we also obtained values (based on radio OH and 
HCN observations and IRAS data) for the total mass loss rate M and the total 
emitted fluxes Fik pc — 47r(l kpc) 2 Ft ota ;; we thus estimated the dust to gas 
ratio fi (in mass fraction) and the distance d for each source. Our results are 
summarized in Table 2. 
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Figure 2: Model spectra for O-rich circumstellar envelopes; the curves have tio 
= 0.001, 0.003, 0.01, 0.03, 0.1, 0.3, 1.0, 3.0 and 10.0; the IR excess increases 
with increasing optical depth. 
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Figure 3: Same as Figure 2, for C-rich circumstellar envelopes having no = 
0.001, 0.003, 0.01, 0.03, 0.1, 0.3, 1.0 and 3.0. 
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-2-10 1 2 

[12]-[25] 

Figure 4: IRAS [25]-[60] vs. [12]-[25] color-color diagram for O-rich (squares) 
and C-rich (triangles) model in Figure 2 and 4 respectively; increasing optical 
depth results in increasing IR excess (from left to right). 
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Note that the dust to gas ratios are generally low (~ 1CT 3 ), with the ex- 
ception of WX Ser. We would like however to stress that estimates of /x are 
affected by uncertainties in the grain average radius a; Eq. 2 shows that Md is 
proportional to a (assuming a correct fitting of Tio), which is not, in this case, 
simply a normalization factor. The distances we obtained are all in accordance 
with estimates by other authors (see e.g. Loup et al. 1993). 

Figure 5 shows the model spectra for a few sources reported in Table 1, super- 
imposed to the IRAS LRS. All sources except those O-rich with less prominent 
silicate features (LML90's "Sil+" and "S" subclasses) can be satisfactorily fitted 
by the available opacities. In the case of the S-star W And ("S" subclass), to 
fit the LRS spectra was necessary to use a mixture of silicate and carbonaceous 
dust in the proportion 1:4, being impossible to model the source using silicates 
dust only. The models for the sources with weak silicate features in subclass 
"Sil+" were obtained by fitting the continuum only; this allows us to produce 
reasonable values for the physical parameters listed in Table 2, but is not able to 
reproduce the correct profile of the dust feature. In our opinion, this is probably 
related to an inadequateness of the available silicate dust opacities for sources 
of LRS subclass "Sil+" and "S" , but we discuss other possibilities in Section 5. 



4.2 Temperature Profiles 

In spherical symmetry, the thermal structure of the model is expressed by a 
radial law T(r) giving the self-consistent dust radiative equilibrium temperature 
for the dust grains. The temperature profiles thus obtained depend on the dust 
chemical composition and on the optical depth. For large envelope radii they 
can be approximated by a power law T(r) cx r~", where a is the asymptotic 
value for the logarithmic derivative of the temperature profile T(r): 

dlogT 

a = hm ~T"j 3 

The inferred a values are restricted to a limited range (see Table 3); however 
those for O-rich sources are systematically lower than those for C-rich ones (0.35 
vs 0.40), while the case of the S-star modelled with a mixed dust composition 
is intermediate. This numerical result is in good agreement with the analitic 
formula derived by Harvey et al. (1991) for spherically symmetric dust shells 
having power-law opacity Q v cx A _/3 and density nir) cx r~ 2 : 

T(r) cx r - 2 /( 4+ « (4) 

For Ossl silicates j3 ~ 2 in the mid- and far-IR wavelength range, while 
for AC1 Rouleau and Martin amorphous carbon is (3 ~ 1 in the same spectral 
range; Eq. 4 would give a ~ 0.33 for silicates and a ~ 0.40 for amorphous 
carbon, similar to our results, appropriate for the outer parts of the envelopes, 
where the mid- and far-IR thermal radiation is emitted. 
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Source IRAS 1 6029-304' 




AFGL '822 







Source IRAS 02143 + 4404 
W Andrornedae 





Figure 5: Model spectra for 6 sources listed in Table 1. IRAS LRS data (trian- 
gles) are superimposed. 
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Name 


Tio 


a 








R LMi 


n i ^ 

\J. ± o 


348 

\J .04:0 


RT Vir 


0.15 


0.349 


S CrB 


0.15 


0.344 


U Her 


0.15 


0.349 


WX Ser 


1.30 


0.347 


AFGL 1822 


8.00 


0.332 


S source 






W And 


0.17 


0.385 


C-rich sources 






Y CVn 


0.01 


0.397 


X CnC 


0.15 


0.400 


FX Ser 


0.20 


0.404 


NSV 11225 


0.20 


0.404 


CW Leo 


1.30 


0.404 



Table 3: Asymptotic values of the thermal structure power-law exponent for the 
modelled sources, ordered for chemical type and increasing optical depth no- 

5 Discussion 

In order to account for the difficulties we found in modelling the O-rich envelopes 
of subclass "S" and "Sil+" , several possibilities can be considered: deviation 
from spherical symmetry, mixing between O-rich and C-rich dust and necessity 
of new opacities for silicates producing LRS other than "Sil" subclass. 

Evidences for non-spherical symmetry of AGB CSEs arise from radio map- 
ping, speckle intcrfcromctry, polarization studies and spectroscopic observations 
as well as from the shapes of planetary nebulae (Balick 1993). Departures from 
spherical symmetry should then be expected also for our sources (see e.g. Silve- 
stro et al. 1996 for a photopolarimetric study of 5 sources in our sample), and 
may affect even the thermal IR emission, as shown in Paper I by the elongated 
shape of CW Leo. Axisymmetric models of AGB circumstellar envelopes have 
been recently developed by Collison & Fix (1991) and Lopez et al. (1995); here 
basically the optical depth depends on the orientation with respect to the line 
of sight. These models thus introduce a new free parameter that allows one to 
fit the dust feature strength; we notice however that the band intensity remains 
unchanged with the slope of the continuum, preventing the correct fitting of LRS 
spectra in our "S" and "Sil+" sources. Furthermore, it is not clear why devi- 
ations from spherical symmetry should be important only for these subclasses 
of envelopes. We conclude that, if a correct two-dimensional treatment of the 
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system geometry is essential for the spatial description of the IR emissions, it 
does not significantly affect the spectra modelling. 

In IE95 the IRAS colors of sources in the 2n class are reproduced admitting 
mixtures of silicates (plus 20% cristalline olivine) and amourphous carbon (in 
the ratio 1:4). The possibility of such a mixing is controversial (dust formation in 
chemical equilibrium models predicts the depletion of the less abundant clement 
between O and C, locked in CO molecules, see e.g. Salpeter 1974), but cannot 
be completely ruled out, since numerous O-rich stars display the SiC feature 
(see e.g. Willems & de Jong, 1986 and Le Van et al., 1992). Even though this 
hypothesis is unattractive, requiring the presence of mixtures only for sources 
of small optical depth, the others being well fitted by ordinary silicate opacities, 
we tested this prescription with our "S" and "Sil+" sources. 

We have not obtained a real improvement in our "Sil+" model with mixed 
dust, because the introduction of a new fitting parameter (the silicate/carbon 
ratio) allowed only to reproduce the shape of the 9.8 /im feature, but not the 
slope of the continuum, even in the LRS region. On the other hand, in the case 
of the "S" source, the injection of a large quantity of amorphous carbon dust 
in the envelope (80% of the total) reduced the strenght of the silicate feature 
making possible to fit the continuum in a large wavelengths range (even thought 
the result of the fitting for the continuum is not as good as in the case of silicates 
only). To understand this result, one should consider that the S-stars like W 
And are considered objects in transition between the O-rich and the C-rich class, 
and thus a contamination of silicate dust in a C-rich circumstellar envelope is 
not at all surprising. 

A third possibility to have weak silicate bands in O-rich envelopes, as re- 
quired in "S" and "Sil+" LRS, without the introduction of carbonaceous dust, 
is to assume annealing and aging of dust grains, as proposed by Nuth & Hecht 
(1990). Their grain condensation scheme starts with oxydation of SiO, AlO and 
OH in O-rich stellar winds, producing grains characterized by weaker features 
and secondary bands at 11-13 /xm; fully grown silicates with strong features 
condensate only later, and are thus present only in more evolved O-rich en- 
velopes. This view is in agreement with the classification schemes of LML90, 
assuming that LRS subclasses are associated with dust of different ages. If the 
Nuth & Hecht (1990) sequence is correct, our O-rich sources with LRS subclass 
other than "Sil" should be associated with envelopes with less evolved dust. 
We have not tested LML90 predictions, because this would have required de- 
tailed dust opacities for each subclass, presently not available. Fully developed 
radiative transfer modelling is in fact necessary even to model optically thin 
envelopes as shown by IE95, and should be performed to extract dust opacities 
from astronomical data. 



15 



6 Conclusions 



Radiative transfer modelling of circumstellar envelopes, based on simple as- 
sumptions as spherical symmetry and steady state outflows, allows detailed 
fitting of infrared data for single sources and the determination of system pa- 
rameters as the distance, luminosity and total amount of dust ejected by the 
central star. The correlation with radio observations (giving informations of the 
total mass loss) will thus provide estimates of the dust to gas mass ratio for 
each modelled source. The accuracy of the fitting procedure depends mainly on 
the available opacities, expecially in the mid-IR range where the dust features 
are present; a good description of the mid-IR spectra is in fact essential for the 
determination of the envelope optical depth and chemical composition. 

The necessity of accurate spectral modelling is becoming more and more 
important taking into account the better spectroscopic capabilities of ground 
based mid-IR cameras and ISO data. 

Our models confirm that the shape of carbonaceous dust features are repro- 
duced assuming a mixture of amorphous carbon, described by AC1 Rouleau & 
Martin (1991) opacity, and (10%, or less) Pegourie (1988) SiC. Minor compo- 
nents (as PAH) might be present, but the quality of IRAS LRS do not allow 
their detection. 

Models for O-rich envelopes associated to M and S stars suggest the pos- 
sibility of different varieties of silicate dust, even if mixing between oxidic and 
carbonaceous dust should be considered, at least in the case of transition objects 
as the S-stars. If the LML90 scheme is valid, we can recognize in the Ossenkopf 1 
"dirty silicates" (OS92) the correct opacity for subclass "Sil" silicates, associ- 
ated to 2n LRS spectra and a mixture of Ossenkopf 1 "dirty silicates" (OS92) 
with a large amount of Rouleau & Martin (1991) AC1 amorphous carbon (in 
the ratio 1:4 for the subclass "S" sources). For the other subclasses accurate 
fits need opacities which at present are not available. 

The dust-to-gas mass ratios we obtained are generally low (between 0.1 and 
0.5 %), except in the case of the O-rich "Sil" source WX Ser (1.5%); there is 
a tendency for having lower dust-to-gas ratios for C-rich envelopes, but this 
result may be influenced by uncertaintes in the average grain radius (that can 
be different between oxidic and carbonaceous dust) and by the inaccurate fitting 
of the silicate features for O-rich sources other than "Sil" subclass. New models 
taking into account the global dynamics of the system and coupling the dusty 
component with the gaseous one will allow direct correlation between IR and 
radio data, and will increase the accuracy in the determination of the dust-to- 
gas ratio, giving better insight on the dust condensation processes responsible 
for the mass loss in AGB stars and the formation of circumstellar envelopes. 
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